Effect of localized holes on the long-range order in bilayer antiferromagnets 
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The effect of localized holes on the long-range antiferro- 
magnetic order in bilayer cuprates is studied, by applying the 
renormalization group to the appropriate classical non-linear 
sigma model. The theory accounts quantitatively for the mag- 
netic phase diagram of Ca doped YBa2Cu30e. 
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The parent compounds of high-T c superconductors, 
La 2 Cu0 4 (LCO) and YBa 2 Cu 3 6 (YBCO), are anti- 
ferromagnetic (AF) insulators. Hole doping suppresses 
rapidly the AF long-range order in both compounds. In 
the monolayer La2Cu04 doped with strontium or with 
excess oxygen, the antiferromagnetic long range order 
disappears at a hole concentration x close to 2% , while 
in the bilayer YBCO the critical concentration, x c , is 
about two times larger, x c « 3.5% [|||- 

In this paper we present the results of a theory, which 
explains quantitatively the strong reduction of the transi- 
tion temperature by doping, and the distinction between 
the magnetic phase diagrams of LCO and YBCO. 

It has been argued in Ref. Q that holes residing on the 
oxygen bonds in the Cu02 planes generate an effective 
Cu-Cu ferromagnetic coupling, which causes canting of 
the spins in the AF background. The canting angle de- 
cays with distance similar to the potential of a magnetic 
dipole. Following Refs. we describe the j'-layered 

system by the reduced Hamiltonian (i.e., the Hamilto- 
nian divided by the temperature T) H = Hq + Hmt- 
Here Ho is the non-linear-cr-model Hamiltonian in the 
renormalized classical region for the undoped sys- 

tem: 



Ho = §; I M£[(VnW) 2 + 7 P (nW-n^ '^ 
J p=i 



) 2 ]}, 



(1) 



where p s is the in-plane spin stiffness, and 7 P represents 
the interlayer coupling, with 71 = 0. 

The Hamiltonian Hi nt describes the interaction of the 
(randomly distributed in both planes) dipole type impu- 
rities with the fields nW @,|§. Denoting by a(rf) the 
unit vector directed along the doped bond at r» p \ and 

by Mm(r < f ^) the corresponding dipole moment (m(r^) 
is a unit vector), we have 



i,p,p' 



(2) 



Here /3 PP > = 5 pp < 
f(p) 



and f> w (r) =M£ ( j(r-r 



i(l - Spp>), T p — 7 P < 1, 
{p) )a t (r ( f ) )m(r ( f ) ). At tern- 
peratures lower than the optical spin-wave gap (approxi- 
mately 70 meV for YBCO M|) only configurations with 

n (l) ~ n (2) 

dominate the action g| , and the Hamiltonian 
transforms to 



(3) 



with t v = T/up s , and ff 3 « 'Epff'/v. 

The randomly distributed effective dipoles develop 
dipole-dipole interactions mediated by the AF spin back- 
ground. Therefore, at sufficiently low temperatures 
they either freeze in a random spin glassy way or at 
least develop exponentially long range spin-glassy cor- 
relations. Assuming that the range of these correlations 
is much larger than that of the AF correlations, which in 
hole doped AF is finite at all T B, justifies treating these 



f(p) 



moments as quenched. The variables r'f and a(rj, p ) are 
also quenched. Denoting quenched averages by [...], we 



write [f^(r)f^(r')]=X u S^6 ij S(r-r'), 
K = M 2 x/6v = A v x. 



with 



(4) 



The ^-layer problem is thus mapped onto that of the 
monolayer one [|]J(|, albeit with rescaled parameters t v 
and A„ . Both these parameters are marginal in the renor- 
malization group (RG) sense. This allows to employ the 
RG methods to find the effect of the dipole moments on 
a two-dimensional (2D) AF. Generalising the results de- 
rived in Ref. || for LCO, with v = 1, the 2D correlation 
length obtained in the one-loop approximation is: 



6d = C{T, x) exp(27r/3A v ), t v < A„, 



(5) 



and 



S.2D = C(T,x)exp(^- 



t y 3£?, 



),A v <t„. (6) 



Here C(T,x) is a smooth function of its arguments, 
known for the bilayer (v = 2) YBCO only at x = 

§• 

The 3D transition temperature T^{x) of a system con- 
sisting of weakly coupled planes may be deduced from the 
relation 



1 



ae 2D (T N ,X)~l, (7) 

with a generated by an interplane (interbilayer) ex- 
change, or some in-plane spin anisotropy. 

Combining Eqs. (||), (||), and (0), we obtain the crit- 
ical line Tn(x). At Tn(x) < p s AiX, the critical line is 
expected to be practically vertical, with x c given by 

x- x = 3A v L(0,x c ), (8) 

where L(T, x) = (l/4?r) \n(a 2 /aC 2 (T, x)), a is the lattice 
constant. At smaller defect concentrations, i.e. forTjv(cc) 
higher than p s Aix, the critical line is 

T N (x) _ A v L(T N ,x)x 

T N (0) ~ 1 - [1 - 3A v xL(T N , x)]V3 ' [ > 



Fig. 1 compares the theoretical critical lines, calcu- 
lated from Eqs. (|p) and (|^) without fitting parameters, 
with the experimental data. The agreement is very good 
for both LCO and YBCO, taking into account the uncer- 
tainties in the parameters A and L, obtained as explained 
above. 

In summary, we applied the renormalization group to 
consider the effect of localized holes (dipole-type impu- 
rities) on the magnetic properties of Ca doped YBCO. 
We calculated the 2D correlation length and obtained the 
phase diagram, in good agreement with the experimental 
findings. 

This project has been supported by the US-Israel Bi- 
national Science Foundation. 
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FIG. 1. Tn(x) /Tjv(O) versus hole concentration x. Full 
lines are theory (Eqs. (8), (9)) for LCO (left line) and YBCO 
(right line). Symbols are from experiments; LCO: • Ref. [12], 
+ Ref. [13] v Ref. [14], A Ref. [3], YBCO: o Ref. [2], A 
Ref. [3]. 
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The parameter Ai sa 20 for LCO was estimated in 
Ref. |(| from the temperature dependence of £ in slightly 
doped samples. It follows then from Eq. (|J) that A 2 for 
YBCO is Ai w 10. Considering that the logarithmic fac- 
tor L does not differ much for different cuprates, we con- 
clude that x c in YBCO should be about two times larger 
than in LCO. Since both TV and x c depend on C{T,x) 
only logarithmically, we neglect in what follows the tem- 
perature and concentration dependence of C (T, x) . The 
prefactor C(0,a; c ) for LCO is approximately 1.26 A [||, 
while a w 10~ 4 §]. This gives L w 0.91. For YBCO we 
estimate L from Eq. (0) for the undoped sample. With 
T N w 410 K 0, and p s k 200 K (this follows from the 
data given in Ref. [O] and the quantum renormalization 
factors given in RefTM), one gets L m 0.97. 
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